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Abstract: An important problem in applied science is the continuous nonlinear filtering 
problem, i.e., the estimation of a Langevin state that is observed indirectly. In this paper, 
it is shown that Euclidean quantum mechanics is closely related to the continuous nonlinear 
filtering problem. The key is the configuration space Feynman path integral representation 
of the fundamental solution of a Fokker-Planck type of equation termed the Yau Equation 
of continuous-continuous filtering. A corollary is the equivalence between nonlinear filtering 
problem and a time-varying Schrodinger equation. 
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1. Introduction 

The problem of the evolution of a Langevin state, or a signal of interest, described by a 
continuous-time stochastic dynamical model arises frequently in practice. Specifically, in 
(classical) filtering probelms one deals with macroscopic objects whose state variables are 
phenomenologically well-described by classical deterministic laws modified by external dis- 
turbances that can be modelled as random noise. So the state of the system is described by 
a noisy version of a deterministic nonlinear dynamical system termed the state model, i.e., 
the dynamics is governed by a system of first-order differential equations in the state variable 
(x(t)) with an additional contribution due to noise that is random {i/{t)). The noise in the 
state model is referred to as the signal noise. If the noise is Gaussian (or more generally 
multiplicative Gaussian) the state process is a Markov process. Since the process is stochas- 
tic, the state process is completely characterized by a probability density function. The 
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Fokker-Planck-Kolmogorov foward equation (FPKfe) describes the evolution of this probabil- 
ity density function (or equivalently, the transition probability density function) and is the 
complete solution of the state evolution problem Q. 

However, in many applications the signal, or state variables, cannot be directly observed. 
Instead, what is measured is a related stochastic process (y(t)) called the measurement pro- 
cess. The measurement process can often be modelled as yet another nonlinear continuous- 
time stochastic dynamical system called the measurement model. Loosely speaking, the 
observations, or measurements, are discrete-time samples drawn from a different system of 
noisy first order differential equations. The noise in the measurement dynamical system is 
referred to as measurement noise. The continuous-continuous "filtering" problem is the es- 
timation of the continuous-time signal or state process given the (noisy) observations of a 
related continuous-time stochastic process (the measurement or observation process). For a 
recent discussion, see ^J. 

The conditional probability density is the complete solution of the filtering problem. 
That is, it embodies all the probabilistic information about the state process contained in the 
observations and in the initial condition. Here the Bayesian point of view is being adopted. 
From the conditional probability density, an optimal estimate may be computed for any loss 
function. For instance, the minimum variance estimate is the conditional mean. The solution 
of the nonlinear filtering problem is said to be universal if the initial distribution can be 
arbitrary^ . 

The traditional approach to the solution of the continuous-continuous universal nonlinear 
filtering problem requires the solution of the Duncan-Mortensen-Zakai (DMZ) equation, a 
stochastic differential equation (SDE) describing the unnormalized conditional probability 
density. The DMZ equation can be gauge transformed to the time- varying partial differential 
equation termed the robust DMZ equation. However, since the robust DMZ equation is a 
partial differential equation (PDF) with the coefficients depending on the measurements the 
PDF cannot be solved off-line, or in real time. In |^], it was proved that solving the robust 
DMZ equation is equivalent to solving a PDF, termed the Yau Fquation (YYe) of continuous- 
continuous filtering, whose coefficients are independent of the measurements. Hence, the YYe 
can be solved off-line and also in a memoryless way. 

Recently, it has been noted that the Feynman path integral^ can be used to solve the 
general nonlinear filtering problem. In fact, it was shown that the path integral formulation 
of the continuous-continuous filtering problem directly leads to the YYe[Q]. The path integral 
formula for the fundamental solution of the YYe was also derived in The advantages 
of the Feynman path integral formulation are many, including a completely independent 
and self-contained formulation and solution of the general nonlinear filtering problem (than is 
traditional in filtering theory literature which is based on measure theoretic techniques termed 

classic discussion of nonlinear filtering theory can be found in the text by JazwinskiQ. For a more 
up-to-date discussion, see 

^For a classic review of path integrals, see [^]. 
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the Feynman-Kac formalism), and simple, efficient and accurate algorithms implementable 
in real-time 

The Feynman path integral has proven to be a very powerful tool in modern theoretical 
physics, often leading to results that are not evident using other methods. In this paper, 
one such instance is provided in the application to the filtering problem. In particular, it is 
demonstrated that there is a very close relationship between nonlinear filtering theory and 
Euclidean quantum mechanics. Specifically, the fundamental solution of the YYe may be 
viewed as the expectation of a certain operator in a Euclidean quantum mechanical system. 
This follows from the path integral formula for the fundamental solution of the YYe derived 
in [Q. This result generalizes the equivalence between nonlinear filtering and the Euclidean 
Schrodinger equation derived in Q via more complicated means. Their result was limited 
to a filtering system where the signal model drift is the same as in the finite-dimensional 
Yau filter^. This result does not require that assumption — the signal model is assumed to 
have additive noise (or somewhat more generally, orthogonal diffusion vielbein so that the 
diffusion matrix is proportional to the identity matrix), and the measurement model noise is 
also additive. Also, no explicit time dependence is assumed in the model. 

The outline of our paper is as follows. In Section ||, notation and some of the important 
results for the real-time solution of the continuous-continuous nonlinear filtering problem 
are reviewed. In Section ^, are summaries of the results obtained using the Feynman path 
integral formulation, specifically the path integral expressions for the fundamental solution for 
the FPKfe and the YYe. In Section ^, the equivalence between Euclidean quantum mechanics 
and nonlinear filtering is presented. In the following section it is shown that it leads to the 
equivalence between nonlinear filtering and the time-varying Schrodinger equation obtained 
in 0. In Section some important conceptual issues are discussed which further clarify 
the relationship between the YYe and Euclidean quantum mechanics. The conclusions and 
directions of future work are presented in Section ^. In Appendix the path integral formula 
for the fundamental solution of the YYe that forms the basis of our work is verified. 



2. Basic Results of Nonlinear Filtering 

A very brief summary of results of filtering theory is provided in this section. For the purposes 
of this paper, the important result is that the solution of the filtering problem requires the 
solution of a Fokker-Planck type of equation called the Yau equation. A more complete 
discussion with references can be found in [^] and It is worth contrasting the simplicity 
of the Feynman path integral solution discussed in the next section with the traditional 
discussion summarized here. 



^Their explicit solution imposed an additional condition on the measurement model drift, but that is not 
relevant for our discussion here. 
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2.1 From the DMZ SDE to the Robust DMZ PDE 



The signal and observation model considered is the following: 

dx{t) = f{x{t))dt + e{x{t))dM{t), x{0)=xo, 
dy{t) =h{x{t))dt + dw{t), y(0) = 0. 

Here x and v are M"— valued stochastic processes, and y and w are R"*— valued stochastic 
processes, respectively, and e G R"^". These are defined in the Ito sense. The processes v and 
w are assumed to be independent Brownian processes with variances hu and respectively. 
Also, f{h) is referred to as the signal (measurement) model drift, e as the diffusion vielbein, 
and ee^ as the diffusion matrix. In this paper, the additive noise case is considered where e is 
the identity matrix, although the same result holds if one assumes that the diffusion vielbein 
is orthogonal. Finally, no explicit time dependence is assumed in the model. 

The unnormalized conditional probability density, a{t, x), of the state given the observa- 
tions {Y{s) : < s < t} satisfies the DMZ stochastic differential equation: 

m 

da{t, x) = ^Y(^{t, x)dt + ^ia{t, x)dyi{t), where (t(0, x) = (To(x). (2-2) 

1=1 

Here 




where ^ is the zero-degree differential operator hi{x),i = 1, . . . ,m, cro{x) is the proba- 
bility density at the initial time to = 0. Under the following gauge transformation 

u{t, x) = exp ^- hi{x)yi{t)^ a{t, x), (2.4) 

the DMZ SDE is transformed into the following time-varying PDE called the robust DMZ 
equation ||8|: 




n(0, x) =ao{x). 

Here A is the Laplacian. The solution of a PDE when the initial condition is a delta function 
is called its fundamental solution. 
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2.2 The Yau Equation of Continuous-Continuous Filtering 

Recently, S-T. Yau and Stephen Yau made a major advance in the real-time solution of the 
general nonlinear filtering problem ||5|. They began by observing that if ui{t,x) satisfies the 
equation 



i=l * i=l \ j=l I 



(2.6) 



dxi 



i=l 



i=l 



i=l i=l 



dhi 



E i;y.(r/)y,(r,)g-(x)|^(x))n,(t,x), 

i,j=l A;=l 

ni(r/_i,2;) =ti/_i(Ti_i, x), 
in the time interval r/_i <t <ti, then the function ui{t,x) defined as 

ui{t,x) = exp ^^yi(Ti)/ii(x)^ ui{t,x) 

satisfies the parabolic partial differential equation termed the Yau Equation (YYe) 
dui 



(2.7) 



dt 



it,x) 



Z^7w^(*'^)-2^/K^)77-(*, 



2 ^ 9x 

1=1 



i=l 



^1=1 * i=i J 



(2.8) 



in the same time interval. The converse of the statement is also true. In j^, they further 
showed that it is sufficient to use the previous observation, i.e., ui{t, x) satisfies Equation 2^ 
if and only if ui{t, x) defined as 



ui{t,x) =exp ^yi(r/_i)/ii(x) ui{t,x) 



(2.C 



satisfies Equation 2^ in the time interval r/_i < t < ti. Equation 2/7 ( Equation |2.9| ) and 
Equation 2.8 are referred to as the post-measurement (pre- measurement) forms of the YYe. 

Equation 2.6 can also be obtained by setting y{t) to y{Ti) in Equation It was proved 
that the solution of Equation 2.f: approximates very well the solution of the robust DMZ 
equation (Equation p.5|), i.e., it converges to u{t,x) in both the pointwise sense and the 



sense. Thus, solving Equation |2.5| is equivalent to solving Equation 2.5. In fact, in |1C] it was 
proved that u{t,x) converges to a{t,x). 

There are two important points to note. The coefficients of the robust DMZ equation 
(Equation |2.5| ) contain the measurements. Consequently, the partial differential equation 
to be solved for continuous-continuous filtering is unknown prior to measurements. In other 
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words, the robust DMZ equation has to be solved on-hne; its solution cannot be pre-computed. 
By contrast, in the YYe (Equation ^I^ ), measurements are absent from the partial differential 
equation. The measurements only enter the initial condition at each measurement step. This 
implies that the YYe can be be solved off-line; there is no need for on-line solution of PDEs^. 
This makes real-time solution feasible even if the state dimension is large. The algorithm 
based on the YYe is termed the Yau algorithm 

In addition, it is noted that all other known algorithms for continuous-continuous filtering 
assume boundedness of the measurement model drift (see discussion in and [0). This 
is a highly restrictive assumption and implies, for instance, that they cannot even handle 
the linear model which the Kalman filter handles very well. The Yau algorithm is shown to 
converge to the true solution under much weaker conditions on the drift (see [^,1^], ||Tl|| ). In 
that sense, the Yau algorithm is the only practical nonlinear filtering algorithm for the general 
filtering problem. Also note that from a path integral perspective, the Yau algorithm is the 
most natural one. 

The simplicity of the YYe has already led to some important advances. In particular, 
it has led to a demonstration of an equivalence between solving a class of nonlinear filtering 
problems and the time- varying Schrodinger equationj^. As a result, the Yau PDE can be 
reduced to a system of ODEs explicitly solvable using the power series method. This result 
will be shown to be a corollary of the main result of this paper. 



3. Path Integral Formula for the Fundamental solution of the YYe 

In recent papers, it was demonstrated that the path integral method leads to an independent 
formulation and solution of the universal nonlinear filtering problem |12, Specifically, 
the path integral formula for the transition probability density in continuous-discrete and 
continuous-continuous filtering was investigated. One of the results of the aforementioned 
work is the path integral formula for the fundamental solution of the FPKfe and the YYe. 
Some of the relevant results from the papers are reviewed below. 

When the diffusion matrix is h,yl, the transition probability satisfies the FPKfe (see, for 
instance, H) 



dp 



(t, x) = - — [fi{x{t),t)pit, x)] + ^Yl 



. , dxi 
1=1 

^p{t,x). 



2 ^dxj 

1=1 ' 



(3.1) 



The fundamental solution of this FPKfe can be shown to be given by (see [13| and []l2| 



P{t,x\tQ,XQ) 



x{t)=x 



x{to)=XQ 



:^x(t)]exp (-— ^ 



dt 



-1 ■^^0 



{±,(t) - fMt))? + K^{x{t)) 



(3.2) 



*This assumes that the measurement steps are equidistant or known. 
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In continuous-continuous filtering, it is necessary to incorporate the measurement stochas- 
tic process as well as the signal process. That is, the measurement process ensemble must also 
be considered. The inclusion of measurement noise means that each system in the ensemble 
leads to a different time-dependent vector y{t). Although only one realization of the measure- 
ment stochastic process is observed, it is still meaningful to talk about an ensemble average of 
the measurement process (in addition to an ensemble average over the state process). Thus, 
the quantity of interest in continuous-continuous filtering is 

P{tiXi;y{ti)\U^i,Xi-i;y{ti^i)) = (^{S'^ixiU) - Xi)5"'{y{ti) - y{U))) |x(t,_i)=x,_i,y(i,_i=y._i), 

(3.3) 

where (•)^ denotes averaging with respect to measurement noise /x(t). It has been shown in 
that the transition probability density conditional on the measurements is given by^ 

P{ti,Xi;y{ti)\ti-i,Xi-i;y{ti-i)) = exp [ ^ hk{x{ti)) [ykik) - ykiU-i)] ) P{ti,Xi\ti-i,Xi^i) 

\^^' k=i J 

(3.4) 



where 



r-x{ti)=Xi / Y \ 

P{ti,Xi\ti.i,Xi.i)= [^x(t)]exp -— 5(i,_i,t,) , (3.5) 



and 



1 



S{ti-i,ti) = - dt 
^ Jti-i 



m 



^(x,,(t) - fMt))? + iSi^it)) + fYl ^i(^w) 

• -, • -, OXi lin 

.«=! 1=1 ^ k=l 



(3.6) 



It can be shown that P{ti,Xi,ti-x,Xi-i) is the fundamental solution of the YYeQ. For 
completeness, a verification of this result is presented in the Appendix. 

The transition probability density is the complete solution to the continuous filtering 
problem. Thus if the initial distribution is u{tQ,x), then the evolved conditional probability 
distribution is given by 

u{t,x)= P{t,x;y{t)\to,x'Q;y{to))u{to,xo){cf'xo} ■ (3.7) 



4. YYe and Euclidean Quantum Mechanics 

In this section it is shown that there is a general equivalence between the YYe and Euclidean 
quantum mechanics. 

^Here the "post-measurement" form is used; the pre-measurement form does not affect the conchisionM]. 
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4.1 A General Equivalence 

Consider the path integral formula for the fundamental solution of the Yau Equation 



P{t,x\to 



i-x{i)=x 
>^o) = / [ 

Jx{to)=XQ 



1 



Mt)]exp{-— 5), 



(4.1) 



where the action S is 

ft 



2 7*0 



dt 



lit ''' M fc 



.fe=i 



j=i 



fe=i 



Integrating by parts, the third term above yields 

" rt 



n „j n px{t) 

J2 dtxi{t)fi{x{t)) = -J2 dxi{t)fi{x{t)) 

i=l •'to j=l •Jx{to) 



(4.2) 



(4.3) 



Here the Feynman convention is used implicitly, i.e., symmetrized arguments of fi, so that 
the ordinary rules of calculus apply. Thus, the action simplifies to 



1 /•* 

S = - dt 

2 7*0 



E m) + fh<t))] ^M^)) + E 



1=1 



1=1 



k=l 



E / dxi{t)h{x{t))- 
(4.4) 



We can therefore write the fundamental solution as 



x{t)=x 



x{to)=X(, 



[^x(t)]exp(-E 



" ^ rx{t)=x 



i=l 



dx- 



x{to)=xo 



:(t)/,(x(i))^exp(^-i-^*dt^j 



(4.5) 



where the "Lagrangian" =Sf , is 

=^=Je-?w+^e 



i=l 



+ |^E^^(-(*))' 



(4.6) 



fe=i 



= T-V, 



and 



r = ^Ex?(t), 



(4.7) 



1=1 

n 



1 



1=1 



5/. 



//(x(i)) + ?i,^(x(t)) 



+ lrE^^(-W)- 



2h 



fe=l 



Consider a Euclidean quantum mechanical system with the Hamiltonian 

je = T + V. 
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Then, performing the path integral quantization of this system wih lead to the following 
expression of the transition probability amplitude: 



(x,t|xo,to) 



x{t)=x 



[&x{t)] exp 



'x{to)=XQ 

Now, let us consider the following matrix element: 



1 



''f J to 



(4.9) 



x{t)=x 



x{to)=Xo 



[^x{t)] exp — 




1 



E 



x{t)=x 



Jx{to)=xo 



dx^{t)fi{x{t)) 



i=l 

xit)=x \ / 

dxi{t)fi{x{t)) exp - — 



x{to)=Xo 



h 



dt^ 



(4.10) 



■u J to 



In the general case, the value of the line integral depends on the path. Therefore, it 
cannot be "factored out". After all, it is a part of the path integral, which is a weighted 
sum over all paths. This is not a "gauge transformation" in the traditional sense (since it is 
not outside of the path integral). This expectation value can be evaluated perturbatively, or 
more generally, non-perturbatively. 

However, it can be viewed as an expectation of the signal model drift line integral oper- 
ator, i.e.. 



exp 



x(t)=x 



x{to)=Xo 



dXi{t)fi{x{t)) 



(4.11) 



in the quantum mechanical system with Lagrangian given by Equation 4.6. 

Note that there are several differences between this operator and the Wilson lines that 
arise in quantum field theories. First of all, this is a quantum mechanical system, not a 
quantum field theory. Secondly, in field theory the coordinates are a parameter, whereas here 
(as in quantum mechanics) the coordinates are quantized and not an operator. 

This can be summarized as follows: The fundamental solution of the Yau Equation, 
Equation 2.8, is the same as the matrix element of the operator in Equation 4-^^ for the 
Euclidean quantum mechanical system with the Hamiltonian as in Equation 

4.2 A Special Case 

The quantum mechanical equivalence can be made even more direct for a certain class of 
the signal model drift. This reault is known and shown here using path integral methods. 
Specifically, suppose the line integral is path independent. Then, it can be factored out of the 
path integral since its value does not depend on the path. It is straightforward to see when 
this is possible. 

In the one-dimensional case, it is always possible to write the (smooth) drift as a gradient: 



(4.12) 



- 9 - 



since g{x) is simply given by the integral of f{x). 

In the general n— dimensional case, path independence of the line integral requires that 
the drift be a gradient: 

/(x) = V0(x). (4.13) 

A special case is the linear symmetric case. Specifically, if L is a symmetric matrix the drift 
as 

n 

f{x) = Y,[LijX,+li], (4.14) 
= V</>(x), 

where 

Hx) = 2 XiLijXj + ^ kxi. (4.15) 

i,j=l i=l 

Clearly, the the state model drift for the Yau filtering system with L a symmetric matrix is 
also a gradient of some scalar. 

In all such cases, the expectation of the relevant operator is simply 

fx(t)=x /if* 



In fact, 



exp (-3- [0(x) - (Pixo)]) I [&x{t)] exp (-^ f dtjA . (4.16) 

/ Jx(to)=xo V J to J 



P(t, x|to, a^o) = exp ( ^ \<^(x) - (^(xq)] ) (x, t|xo, to) • (4.17) 



A gauge transformation relates the fundamental solution of the Yau Equation and the corre- 
sponding Euclidean quantum mechanical system. 

The discussion in this subsection can be summarized as follows: When the signal model 
drift is a gradient of a scalar field, the fundamental solution of the Yau Equation is, up to a 
gauge transformation, the same as the transition probability amplitude of a Euclidean quantum 
mechanical system with Lagrangian given by Equation |^.6|. 



5. Nonlinear Yau Filtering System and Time-Dependent Schrodinger equa- 
tion 

5.1 The Work of S-T. Yau and Stephen Yau 

As discussed in Section |2[ in order to obtain the unnormalized conditional probability a{t, x), 
it suffices to solve the YYe 

^{t,x) = lAn(t,x) - f;/,(x)^(t,x) - ^M(x) + if;/i2(x) n(t,x). 



j=i "J \i=i ^ j=i 

u(ri_i,x) = ai{x). 
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In the authors considered the filtering system signal model drift 

f{x) = Lx + l + V<l), (5.2) 

where L is an anti-symmetric matrix. The symmetric part can be incorporated in (j). By 
considering the quantity ^{t, x) defined by 

=e'^(^)p(t,x), (5.3) 
= e'^(^)i/(t, x), i{t) = B{t)x{t) + 6(t), 

they showed that it is sufficient to solve the following Schrodinger equation 

ff (t,x) = \Vv{t,x) - \q {B-\t)x - B-\t)b{t)) u{t,x), 
U{Ti-i,x) =a,(x)e-*(-), 



(5.4) 



where 



and 



Here, 



q{x) = V20(x) + |V(/>p(x) + 2{Lx + /) • V(/> + ^ hj{x) + tr L, (5.5) 



i=l 



B{t) = e-^\ b{t) = - [ e-^nds. (5.6) 

Jo 



^{t) = -B{t)L, and ^^(t) = -B{t)l, (5.7) 
and B{t) is an orthogonal matrix. 
5.2 Path Integral Derivation 

We now derive using path integral methods the result discussed in the previous sectoion 
relating the nonlinear Yau filtering system to a Schrodinger equation using path integral 
methods. In matrix notation, the Lagrangian for the Yau filter state model is 



1 h h 

^=-{x-Lx-l-\74>)' + ^\/-f + ^^hlix). (5.8) 

^ k=l 

Since 

V • / = tr L + VV, (5.9) 



it follows that 



11 h 
^ = -{x-Lx-lf + -{V(t>?{x) -x-V4> + {Lx + l)-V(f> + ti L + S7^cj) + -^y] hl{x) 

(5.10) 



1 2 1 

-{x-Lx-l) - X ■ 'V4> + -q{x) 



-11- 



where q{x) is as defined in Equation |5.5| . The x ■ V4>{x) term can be integrated out of the 
path integral and simply yields a phase factor 

exp - cPixo)]] . (5.11) 



''U 



Now, in terms of x{t) 



xit) = B-^{t)[x{t) -b{t)], (5.12) 
x{t) = B-^{t)[3:{t) + Li{t) + B{t)l - Lb{t)], (5.13) 
Lx{t) = B-^{t)[Lx{t) - Lh{t)] 



so that 



x{t) - Lx{t) - I = B-^{t)i{t), and (5.14) 
(x-Lx-/)2 = 5^ (5.15) 
Therefore, in terms of x{t), the Lagrangian is 

^=h' + lq{B~Ht)[S:-b{t)]). (5.16) 
Combining this with Equation 5.11 , the fundamental solution becomes 

P{t,x\to,xo) = e^p(^[^{x) -(Pixo)]) [ ''^ [&x{t)]e^p(-^ [ dtjA, (5.17) 
yriu J Jxiito)=xo \ "''^ J to J 

= exp {^[(t^ix) - (t){xo)]^ i'{t,x\tQ,XQ), 

where is given by Equation |5.16| , and this equation defines i>(t, x|to, a^o)- Note that the 
Jacobian of the transformation from the measure [^x(t)] to [^x(t)] is unity since B{t) is an 
orthogonal matrix. 

It follows from the standard path integral formulation of quantum mechanics that v{t,x\tQ^ xq) 
is just the path integral representation of the fundamental solution of the following Schrodinger 
equation: 

^(i, x) = -VMt, - -q {B-\t)i - B-\t)b{t)) u{t, x). (5.18) 
This is precisely the equation obtained in Q. For the initial condition, since 

u{t,x)= / P{t,x\to,xo)a{to,xo)ctxo, (5.19) 



exp ( ^'/'(a;) ) i>{t,x) 



it follows that 



exp ^^0(x)^ j i'{t,x\to,xo)e-K.p ^-^(t){xo)^ a{to, xo)(r-xo, 



P{t,x\to, xo) = exp ( ^ [(p{x) - (pixo)] ] v{t, x|to, a^o), (5.20) 



which is the same as Equation |5.17| . Thus the result in Q has been independently established. 
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6. Additional Remarks 



A few comments on our results are in order. 



1. The general equivalence proposed in Section 4J is that between the fundamental so- 
lution of the FPKfe and a certain matrix element of a quantum mechanical system. 
Note that all previous equivalences were between fundamental solutions of FPKfe and 
a Schrodinger equation. While the previous equivalences were obtainable using opera- 
tor methods, it is not clear how the proposed general equivalence can be derived using 
operator methods. 



2. The Feynman path integral methods used here are very different from the measure- 
theoretic methods used to study the nonlinear filtering problem (see, for instance, |14|). 



The Feynman path integral approach developed in |Q, |2j and this paper leads to an 
independent way of tackling the universal nonlinear filtering problem. In particular, 
unlike the standard filtering theory approaches, the DMZ equation (or its variants) is 
not taken as an input. On the contrary, the YYe is obtained directly as a consequence 
of the path integral approach. It is also noted that the Euclidean quantum physics 
referred to in the equivalence to nonlinear filtering developed in this paper is a quantum 
mechanical one, not a quantum field theoretical one. Also note that hi, here is analogous 
to the Planck's constant, h, in quantum physics. 



3. Although a formal equivalence has been developed between Euclidean quantum me- 
chanics and universal nonlinear filtering, it is important to point out that there is a 
profound difference between classical and quantum probabilities. In the "real time" (as 
opposed to Euclidean time) quantum mechanics, the transition probability amplitude is 
not a probability; it may not even be real. The probability amplitude is to be multiplied 
with its complex conjugate to obtain a probability density. In contrast, the transition 
probability density in filtering theory is a classical probability density. We have shown 
that, in a mathematical sense, matrix elements in a Euclidean quantum mechanical 
system equal the transition probability density of a classical stochastic process. This 
equivalence is purely mathematical, not conceptual. 

4. The reason for some of the mathematical equivalence between nonlinear filtering and 
quantum physics is the semi-group property. That is, in stochastic processes the 
Chapman-Kolmogorov semi-group property is a fundamental property of Markov pro- 
cesses. The Chapman-Kolmogorov semi-group property allows the transition probabil- 
ity density in stochastic process to be written as follows. Let us partition the time 
interval [to,t] into equi-spaced time intervals so that tj = to + ie where e = {t — t())/N. 
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Then, from the Chapman-Kolmogorov semi-group property it follows that 



■x{t)=x 



P{t,x\to,Xo) 



Kx(ti) • • • d''x{tN-l)}P{t, x\tN-Ux{tN-l)) ■ ■ ■ P{ti,x{ti)\to, Xo) 



x(to)=Xo 



(6.1) 



„x{t)=x (N~l "I N 

Jx{to)=Xo [ J 




where x(to) = and x{t]\f) = x and the delta function condition in the definition of 
P{t, x\tQ,xo) is written as integration limits. On the other hand, in quantum physics, the 
semi-group property is a basic property of the time-evolution operator. It is the semi- 
group property of the time-evolution operator that leads to the path integral formula 
for the transition probability amplitude. 

5. Note that the unitarity of the time evolution operator plays a crucial role in quantum 
mechanics. This is because it is essential for conservation of probability. It is noted 
that the unitarity of the time evolution operator implies that the Hamiltonian operator 
is Hermitian®. 

The evolution of the probability distribution for a continuous Markov processes is de- 
scribed by the FPKfe. The FPKfe operator is not a Hermitian operator. This is not 
inconsistent with the conservation of probability. This is because the FPKfe is a con- 
tinuity equation, and so the probability is conserved (as long as the boundary terms 
vanish). This is yet another instance of the profound difference between classical and 
quantum probabilities. 

6. In the continuous-continuous filtering case, the YYe plays the role that the FPKfe plays 
in continuous-discrete filtering. However, the YYe is not a continuity equation. This is 
not a contradiction since the YYe is related to the unnormalized conditional probability 
density, whereas the FPKfe evolves (and preserves the normalization of) the normalized 
probability density. 

7. In this paper, it has been assumed that the noise is additive and the model is not 
explicitly time-dependent. In the more general case, a simple equivalence, as derived in 
this paper, is not possible. For instance, in order to obtain the YYe, it was necessary 
to assume that the measurement model was not explicitly time dependent (see Q ) ; this 
is not valid in the general case. Also, when the noise is multiplicative, quantization 
is not as straightforward. This can be traced to the well-known operator ordering 

''In the standard formulation of quantum mechanics, Hermiticity of the Hamiltonian is required in order 
to ensure that the eigenvalues of the Hamiltonian (and hence the possible energies) are real and that the 
time evolution is unitary (i.e., conserves probability). This can also be ensured even if the Hamiltonian is 
not Hermitian provided that the Hamiltonian is space-time (or ^^) reflection symmetric. For a pedagogical 
discussion of non- Hermitian quantum mechanics, see [[L5[. 
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ambiguity in quantum physics (since the position and momentum operators do not 
commute). Furthermore, standard calculus manipulations are no longer possible in the 
path integral. However, the path integral result is the same even for the multiplicative 
noise case if the diffusion matrix is proportional to the identity matrix. 

Finally, note that the results of this paper also apply to the FPKfe itself, which is simply 
the case h{x) = 0, in a straightforward manner. 

7. Conclusions and Future Work 

The main conclusion of this paper is that certain continous-continuous nonlinear filtering 
problems are related to Euclidean quantum mechanics. Specifically, the transition probability 
density for nonlinear filtering problems with additive noise and with square diffusion vielbein 
and explicitly time-independent drift is the fundamental solution of the YYe, and is the same 
as the expectation of a certain operator in a quantum mechanical system. A corollary of this 
general result is a derivation of the relationship between the YYe in the Yau filtering system 
and a time-varying Schrodinger equation which was first derived earlier. This equivalence 
leads to some useful results from the methods used in quantum physics. 
There are many possible directions for future work: 

1. The path integral formula can be exactly solvable, or lead to simpler methods of ob- 
taining solutions, in some simple cases. In fact, for the model studied in 0, an explicit 



path integral solution is described in a paper currently in preparation |1£] 



2. The equivalence to Euclidean quantum mechanics immediately leads to many filtering 
problems that can be solved exactly, as will be explored in future papers. The re- 
markable fact is that many exactly solvable Euclidean quantum mechanical problems 
correspond to filtering problems that are not finite-dimensional (i.e.. Lie algebra of 
and hi{x),i = 1, . . . ,m is not finite-dimensional, see |^) . Thus, simplicity in filtering 
theory does not imply finite dimensionality. 

3. The path integral formulation naturally leads to a perturbative solution of the nonlinear 
filtering problem. Such a perturbative solution which is analogous to extended Kalman 
filtering is called extended Yau filtering (EYF). Thus, it is possible to perturb about the 
Yau filter (a generalization of the linear Kalman filter), rather than the linear Kalman 
filter. Clearly, since the EKF is a special case of the EYF, such a formulation will be 
superior to the EKF. 

4. However, it is noted that the fundamental solution is defined nonperturbatively. This 
is important since sometimes the perturbative approaches, like EKF, fail. 

5. Perhaps the most important advantage of the path integral lies in numerical methods of 



computation [17|. In fact, path integrals are the only known way to carry out nonper- 



turbative computation in quantum chromodynamics (QCD). Note that QCD, the gauge 
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theory of strong interactions, is a quantum field theory, not merely quantum mechanical 
and the QCD Lagrangian is highly nonlinear. The excellent numerical results in QCD 
suggest that currently known path integral methods should give very good performance 
for the simpler case of (large dimensional) universal nonlinear filtering. It is sufficient to 
note that the Dirac-Feynman approximation, the crudest approximation of the path in- 
tegral, already yields excellent results (see |Q] ) , and is adequate for smaller dimensional 
problems. 



It is clear that many lines of investigation are suggested by the path integral methods 
and it is planned that results of those investigations will be presented in future papers. 



A. Verification of the Path Integral Formula 

It is now proven that the path integral formula Equation p.4| satisfies the YYe. This closely 
follows the method used by Feynman to verify the path integral formula for the Schrodinger 
equation in quantum mechanics. 

According to the Chapman-Kolmogorov semi-group property 

P{t + e, x\tQ, xo) = j P{i + x')P{t, x'|io, a;o){d"a;'}, (A.l) 

where P{t,x'\tQ,XQ) is given by Equation |3.5| . When e is an infinitesimal 

P{t + e,x\t,x')= (A.2) 
— \xi — x\ — efj(x)l^ — e- '-J-^{x) — e | , 

4=1 «=1 ^ k=\ / 



exp 



2 

satisfied: 



where x = ^(x + x'). The dominant contribution occurs when the following condition is 

x-x' - ef{x) « 0. (A.3) 



We may write in this region 



x = x' + ef{x) + rj, or (A.4) 
X = x' + e/(x) + r/. 



where the equalities are valid to 0(e). Substituting this into Equation A.2| (keeping terms up 
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to 0(e)) 

Pit + €,x\to,xo) = A[ exp(--^^r/2| (A.5) 



2Ke . ^ 

1=1 



m 



i=l * ^ k=l / 



A I exp 



- e/(x,t) - 7]\to,xo). 

It is noted that the Jacobian of the the transformation from xq to rj to 0(e) is included in 
the second step. 

The constant A is fixed by 



A 

Hence, it fohows that 



The left hand side of Equation A.5 



i=l 

is 



dP 

P{t, x\to,xo) + €—{t,x\to, xo). (A.8) 
The second order expansion (in ij) of P{t, x'\tQ, xq) yields 

P{t, x'\to,xo) = P{t, X - ef{x, t) - T]\to, xo), (A.9) 

(" 9 f 1 ™ \ 

1-^Eaf (^'*)-^^E^fc(^) 
i=l * ^ k=l / 

" dP 1 " d^P 

We are interested only in terms of 0(e). The only nonvanishing terms are terms linear in e 
and quadratic in r]. Then it is easy to verify that the diffusion part of the YYe follows from 
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the term quadratic in i] and the rest fohows from the term linear in e, so that 



dP 
'dt 



{t,x\to,xo) 



h d'^P " d r 



2 ^ dxl 

1=1 ' 



fi{x)P{t,x\to,xo) 



i=l 



P{tQ,x\to,xo) = 5"-{x - xo). 



Hence, P satisfies the YYe with a delta-function initial condition. 
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